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Abstract. We obtain a well-behaved continuum limit of projected entangled pair states 
(PEPS) that provides an abstract class of quantum field states with natural symmetries. 
Making use of the recently introduced path integral representation of one-dimensional con- 
tinuous matrix product states (cMPS) for quantum fields, we demonstrate how symmetries 
of the physical field state are encoded within the dynamics of an auxiliary field system of 
one dimension less. In particular, the imposition of euclidean symmetries on the physical 
system requires that the auxiliary system involved in the class' definition must be Lorentz 
invariant. The physical field states automatically inherit entropy area laws from the PEPS 
class, and are fully described by the dissipativc dynamics of the lower dimensional auxiliary 
field system. 



1. Introduction 

The quantum states we observe in nature are highly atypical as compared to a state ran- 
domly chosen from the full configuration Hilbert space % [16]. Indeed, observable states 
comprise only a tiny submanifold of Ti — the physical corner of Hilbert space — whose 
points exhibit highly nongeneric features such as nontrivial clustering of correlations and en- 
tropy areas laws 13] . It is extremely desirable to develop an efficient parametrisation of this 
manifold as this would considerably ameliorate the computational costs of solving physical 
models and provide new analytical tools for the study of quantum field systems. Indeed, 
even a partial parametrisation of the physical corner provides a powerful tool as it supplies a 
variational class useful for the description of low-energy physics. 

The canonical example of such a class of quantum states appears in the setting of one- 
dimensional lattices. There matrix product states (MPS) [I] have enjoyed remarkable success, 
not simply for the calculation of physical properties of strongly interacting lattices, but also 
for such things as the classification of quantum phases, providing a natural foliation of states 
in terms of entanglement, and the construction of exactly solvable models [23l |27J [21] . It 
is also well-established that MPS satisfy two important criteria. Firstly they constitute a 
complete class of quantum states, in the sense that by increasing a "bond dimension" D one 
can capture any pure quantum state of the system. Secondly the class is efficient in the 
sense that the computational cost of calculating expectation values scales polynomially in 
the number of variational parameters. 

The MPS class has provided a fruitful basis for generalisations: by understanding the 
structure of quantum entanglement in such states they have inspired several powerful exten- 
sions to higher dimensions and different geometries. The two most prominent examples for 
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higher-dimensional lattice systems are the projected entangled-pair states (PEPS) (T9| and 
the multiscale entanglement renormalisation ansatz (MERA) [2U [22] • Both of these varia- 
tional classes have proved invaluable in the investigation of strongly correlated physics. So 
far, however, all of these results have been restricted to the lattice setting; the study of con- 
tinuous quantum systems using these classes has traditionally proceeded by first discretising 
the system on the lattice and then employing them as a variational ansatz. 

Continuum systems bring considerable difficulties when it comes to variational computa- 
tions because optimisations can be dominated by UV physics at the expense of infra-red 
physics which ruins the estimation of observables of physical interest [5] . Both of these diffi- 
culties have been addressed with the introduction of special continuum versions of the MPS 
and MERA classes [221 HH IB]- The cMPS class is remarkable in that it requires (in the 
translation-invariant case) only a finite number of variational parameters to specify, but is 
expected, by analogy with the discrete case, to be both efficient and complete in the sense 
already described. Further studies also indicate that cMPS and cMERA are not disrupted 
by the presence of UV divergences [H IB] • 

The MPS and cMPS classes are restricted to one spatial dimension. Therefore it is natural 
to contemplate generalisations of the PEPS class to the continuum. A first step in this 
direction was taken in [14], where such a generalisation was proposed. Unfortunately the 
formulation presented there does not make symmetries manifest and, further, has no obvious 
relationship to the PEPS class. 

In this paper we make use of a recent representation of one-dimensional cMPS states in 
terms of path integrals [2] as a guide for the generalisation of the class to higher dimensions. 
Specifically, we start with a particular sequence of lattice PEPS with decreasing lattice spac- 
ing and obtain a continuum limit using a path-integral representation. This generalisation 
takes the same functional form as the one- dimensional cMPS and manifestly exhibits spatial 
(i.e., rotational) symmetries. The derivation via a sequence of lattice PEPS means that the 
resultant class of field states automatically obey entropy area laws. Furthermore, the impo- 
sition of spatial symmetries on the physical field state is obtained by encoding the symmetry 
into the dynamics of an auxiliary boundary system with the novel result that the dynamics of 
the boundary system is given by the imaginary time evolution of a Lorentz invariant system 
of one lower spatial dimension. We build directly on [2] and therefore assume the reader is 
familiar with the contents of that paper. 



2. Path integrals for cMPS 

Continuous matrix product states are a variational class of states for one-dimensional 
quantum fields [2"2"l [T_3] U\- The simplest formulation is for the single-species bosonic case 
on the real line with field annihilation and creation operators ^a{ x ) an d wa( x ) such that 
[iPa{ x ): VuG/)] = $( x ~ v)- A cMPS is then defined in terms of the quantum field A and a 
(bosonic) auxiliary Z)- level quantum system B by 



|X> = (uj L \Vexp 


r'J 


f K(s) (8)1 + iR(s) ® $ A (s) - iR j (s) <g> ^(s) ds 
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where K is a D x D hermitian matrix and R is D x D complex matrix, \ool,r) ar e D- 
dimensional states of the auxiliary system B, ip A {s) is a bosonic field operator on the physical 
system A, \fi A ) is the Fock vacuum, and V denotes path ordering. 

We now recall the main result of [2] : we can represent a general cMPS via the path integral 
over an auxiliary (0 + l)-dimensional D component complex field 



(2.2) 



V 2 <J) exp [iSfatf)] |$), 



where the path integral is a coherent-state path integral over all configurations of the D- 
dimensional complex vector 0, |$) is a physical field coherent state 



(2.3) |$) = exp 
$(s) = (j)'R(j), and the complex action S is given by 

(2.4) S(4>,ft) = f ds(i 



^*(s)tp A (s) ds 



Kef). 



The case where the auxiliary system is fermionic is identical, except that is a vector of 
Grassmann fields. The expression (12.21) makes equal sense for quantum systems on [0, 1] as 
for the infinite case (— oo, oo). 

The path integral representation of |%) is extremely suggestive of a higher dimensional 
generalisation, namely, we should have that 



(2.5) 



\ X )= / P 2 0exp[-5(0,0t)] |$), 



where the path integral is now over an auxiliary (d — 1 + l)-dimensional field 
components, |$) is now a higher-dimensional field coherent state 



(z, t) with D 



(2.6) 



|$) = exp 



$(x)^(x)-$*(x)^u(x) dx 



\n A ), 



whertj] $(x) = ft(x)R(j)(x), R is a D x D matrix, and S is a local action for a D component 
auxiliary boundary field <fi living on an auxiliary boundary space of one lower dimension. 
Note the replacement of the i with —1 in the exponential of the action: the euclideanised 
formulation will prove vital when imposing symmetries on the physical state. 

3. The continuum limit of a PEPS class 

While the continuum limit of the one-dimensional MPS class is comparatively straightfor- 
ward, the two-dimensional equivalent has proved more challenging. It is true that we can 
simply posit the form of a two-dimensional (or higher-dimensional) cMPS as being generated 
by the continuous measurement process of a lower-dimensional auxiliary boundary field [T4] , 
however this is unsatisfactory for at least two reasons. Firstly, in such a setting it is not clear, 



^Here the notation 0(x) means that we evaluate (j> at the spacetime location (z, t) found by regarding the 
first d — 1 components of x as spatial coordinates and the dth component as a temporal coordinate, i.e., 



3 =0,1, 



, d — 2, and t = x&-\. 
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a priori, how one might impose certain desirable symmetries, such as rotational symmetry, 
on the physical quantum state. Any variational class intended for the efficient description of 
real-world physics should be capable of manifestly exhibiting such symmetries. Secondly, for 
discrete systems higher- dimensional generalizations of MPS already exist, such as the PEPS 
class, which have been powerful tools in understanding the physics of local hamiltonians. As 
such it is also of theoretical importance that we arrive at a continuum limit of PEPS that 
mirrors the one- dimensional cMPS class. 

In [2] a path integral representation was obtained for the one-dimensional cMPS class 
from the traditional discrete MPS class by taking a well-behaved continuum limit. It is this 
representation of a cMPS that we use as our guide for constructing higher- dimensional classes 
with manifest symmetries; in field theory, path integral formulations are ideally suited for the 
imposition of symmetries that would not be manifest according to, e.g., canonical quantisation 
of the field. Our strategy is then to develop a continuum limit of of a special class of PEPS 
such that the physical state is a superposition of field coherent states with amplitudes given 
by a path integral over an auxiliary system and such that the desired symmetries are manifest. 

3.1. The basic Tensor Network setting. A natural higher- dimensional generalization of 
MPS are the Projected Entangled Pair states (PEPS), which are examples of Tensor Networks 
[23| IT2| [20] . The original formulation of PEPS rested on distributing maximally entangled 
pairs of D-dimensional quantum systems between neighbouring sites on a graph, and then 
locally mapping the systems at each point into a single d- dimensional Hilbert space. The 
PEPS construction for arbitrary D can describe any quantum state, and is naturally suited 
to systems displaying area laws. A generic PEPS has an expansion in terms of a product 
basis with expansion coefficients given by a contraction of tensors with respect to a 

particular graph T(V, E): 

(3-i) ix> = c^U)--^--- A ;U)]^-'-^- 

where C denotes a complete contraction of the auxiliary indices (i ■ ■ ■ k) according to the graph 
edge structure E, and r±, . . .r^ label the (product) configurations of a the discrete physical 
systems located at each vertex of the graph. 

An initial instinct would be to begin with a two-dimensional square lattice, and embed 
the discrete system into the one-particle sector of a system of bosonic or fermionic auxiliary 
fields, as was done previously for the 1-D cMPS path integral. If one directly follows this path, 
passing from the discrete PEPS to a continuum path integral, one finds that the underlying 
square lattice structure persists in the field, and one does not obtain a rotation invariant 
physical state (see Appendix EJ. Here we adopt a slightly more involved strategy to handle 
this unwanted feature. 

We begin with a graph of N 2 points {(x,y) = (ne,me) : n,m integers}, with a physical 
spacing e. We view the y spatial direction as an auxiliary time t. As written (13. ip involves 
a contraction over N 2 auxiliary subsystems distributed over the graph, with independent 
couplings to the physical degrees of freedom at each site. Since we wish to view the y 
direction as an auxiliary time dimension we regard the contraction over the N 2 subsystems 
as the sum over configurations of N auxiliary subsystems subject to a sequence of N dynamic 
transformations. The upshot is to replace the contraction along the y direction with a product 
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of N square matrices of dimension D . This is well-known as a transfer operator approach, 
i.e., we simply view the contraction from one value of y to the next as multiplication by a 
transfer operator. The state (13. ip can then be written as 

\x) = sMu(t = i 7 t = o)\u R ) B \n) A , 

where U is an operator on AB given by the (time-ordered) product of transfer operators 
local with respect to the graph sites and \fl) A is some initial product state of A. Specifically, 
U = TU l t=o[M(t)} = M(l)M(l - e)M(l - 2e) • • • M(0), where the transfer operator M(t) 
generates an elementary time-step of size e and is built from local operators on A and B. We 
then follow the idea of [2] and describe the auxiliary systems at each of the lattice sites of B 
with the single-particle sector of the Fock space built from C D . 

The contraction of indices depends on the particular graph structure being useojE However, 
our goal is to construct cMPS states with symmetries and we follow the key principle that 
the symmetries of the physical state are encoded in the dynamics of the auxiliary system. 
For example, a natural symmetry to demand is that of rotation invariance in the spatial 
coordinates of the physical field state. Assuming a state of the form ( 12. 5 p implies that the 
auxiliary action S is invariant under the induced SO(2) rotation group (assuming the measure 
is also invariant). By demanding that the auxiliary system is a physical system we deduce 
that S should be an action describing the completely positive dynamics of an auxiliary system 
(after we trace out A). However, encoding the symmetry into the dynamics of the auxiliary 
system means imposing invariance under SO(2). This implies that the dynamics should 
be viewed as the imaginary-time evolution of a Lorente-invariant system (which is still a 
completely positive map of the quantum state). 

Our strategy is then to first construct a Lorentz-invariant auxiliary action from the con- 
tinuum limit of a sequence of discrete PEPS. We then construct an analytic continuation 
to the euclidean setting and obtain a one-parameter family of discrete PEPS states giving a 
rotation invariant euclidean action as e — > 0. There are clearly different possible choices for 
a Lorentz-invariant action; motivated by the first-order action ( 12 .4p . and certain convenient 
properties of coherent field states, we derive a Dirac-like action from a specific sequence of 

2 It is useful to emphasize that technical subtleties arise when taking the limit of lattice systems. Specifically, 
we might consider a family of graphs {A/c} indexed by some variable k = 0, 1, 2, ... , that converges to some 
dense subset of a compact spatial region A. To each point x £ A& we have an associated Hilbert space 'H(x), 
which could be a space of finite or countably infinite dimension. The total Hilbert space for the full graph 
system is then given by Jik = <8> x eA fc ?^( x )j and in the thermodynamic limit k — > oo, the resultant space will 
have an uncountable dimension. It is common to instead work with a much smaller, separable Fock space 
3t%) constructed to ensure that every state in ^(H) has finite particle expectation value, and splits up into 
a sum $(H) = ® n n {n) of particle sectors H {n) with finite particle numbers. Central to the formation of 
this Hilbert space is the identification of a vacuum state, from which the different n-particle spaces H^ n ' arc 
obtained through the action of creation operators obeying the desired statistics. It is well-known that the 
Stonc-von Neumann theorem fails for these systems, and many unitarily inequivalent Fock spaces may be 
constructed through the choice of different vacua and creation/annihilation operators. For our analysis of the 
discrete to continuum limit, we specify the local Hilbert spaces at each point on the graph, but ultimately 
we make use of a Fock space construction for the state \x) £ 'S('Ha)> an d work with a particular choice of 
creation/annihilation operators for both A and B, with the auxiliary system B carrying bosonic or fermionic 
statistics. 
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FIGURE 1. Tensor network structure for the transfer operator generating one 
temporal layer of the PEPS sequence: here % a b denotes the tensor product 
"Ha <8> Hb for the internal degrees of freedom of the auxiliary system H A . 



PEPS. (We will discuss general euclidean actions, including second-order actions for bosonic 
auxiliary fields, in a future paper.) 

Our first task is to arrive at a Lorentz-invariant situation, a problem for which physical 
intuition is readily available. Since we are looking for an auxiliary (1 + l)-dimensional lattice 
system with locally defined dynamics we assume that each site (x, t) has contraction links to 
future sites (x,t + e),(x — e,t + e) and (x + e, t + e) and also to past sites (x, t — e) , (x — e, t — e) 
and (x + e, t — e). The simplest such choice is to build the operator M out of quadratic terms 
involving creation and annihilation operators; to arrive at states with rotational symmetries, 
we can also make use of spinorial expressions. To generate the spinorial structure we assume 
that at each site x, in addition to the 'flavour' indices i, j,k, . . . , we have access to two internal 
degrees of freedom, with annihilation operators ak, x and bk )X at each spacetime point. We 
also note that the bosonic and fermionic cases can be treated simultaneously by being careful 

with the ordering of terms. Thus we have [a,j iX , a>k,y]± — [ a ],xi a I,J± = 0' i a j,x, a k y ]± — 8jk8 x ,y, 
with similar expressions for b, where ± labels the choice of bosonic or fermionic auxiliary 
system. 

3.2. Generation of a kinetic term and flavor doubling. We now construct a transfer 
operator M e = 1 + eH tot from terms quadratic in a and b, and generate a PEPS via 

(3.2) \ X e) = B (wL\M e (!)M e (! - e)M e (l - 2e) • ■ ■ M e (0)\u R ) B \Q) A , 

in such a manner that the continuum limit has the desired Lorentz symmetry. It should 
perhaps be emphasized that the generators that we construct relate entirely to the auxiliary 
system A, which only forms an abstract device to describe the physical state |% e ), and does 
not commit us to any particular realisation for the two-dimensional physical field system B. 
However, by treating the auxiliary system as physical we can make use of natural intuitions 
of particle interactions when we construct the abstract PEPS class through i^tot- Certain 
assumptions are natural to impose on the terms appearing in H tot , such as left-right symme- 
try, symmetry between a-particles and b-particles and conservation of total particle number, 
however the key term in the construction is the a nearest neighbour 'hopping' term, which 
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we take to be 

(3.3) H h (t) = -^2 J jk (t)( a i x b k , x -e + a\ x b Kx + at x b kjX+e + h.c), 

X 

where x runs over N spatial points, and J J (t) measures the strength of the spatial hopping, 
which for simplicity we take as constant along the spatial direction. 

To analyse this, we perform a discrete Fourier transform in the spatial direction to obtain 

(3.4) H h = - J jk {t){l + 2 cos pe){a\ p b k , p + b] >p a k , p ), 

where a JiP = ^ x e~ tpx a,j jX is the Fourier-transformed annihilation operator (similarly for 

bj iP ), and p = 2im/Ne runs over N points in the reciprocal lattice, for n an integer. 

Our concern is the continuum limit, e — > 0, where the dominant contributions of Hh will 
come from the "low-energy regime" of momenta p near to the zeroes of (1 + 2cospe). These 
occur at the points g M = (— l) M (27r/3e), for /i — 0, 1. The contributions from the two points 
give rise to two flavors in a similar manner to fermionic doubling on the lattice, however it 
must be noted that this doubling only occurs for the auxiliary system, and so the physical 
system (whether bosonic or fermionic) is unaffected. Overall, in the continuum limit Hh splits 
into Hh = Hhfl + Hh ; i with the contributions from the two decoupled flavors given by 

(3-5) HhAt) = V3 (-lYJ j \t) b\ qii+p )pa, 

bl<l/e 

We can redefine J jfc — > J jfc / \/3 and relabel the mode operators as Q>(j,n);p — ®j,qu.+p an d also 
b(j,n); P = bj,q^+p to obtain 

( 3 - 6 ) H h®= E Jjk ^ H»,p \,),v)v<^< 

|p|«l/e;/4=0,l 

This term alone would produce ID massless Lorentz-invariant flavors labeled by the com- 
pound index (J,fj), in which the (J, 0) -field is related to the (j, l)-field through the discrete 
parity transformation V, given in 1+1 dimensions as V = 7 = a z , which inverts chirality. 
This is consistent with the Nielsen- Ninomiya theorem [13], which requires doubling in order 
to achieve a translationally invariant spinor action with chiral symmetry in the continuum 
limit of a lattice model. 

3.3. Coupling and decoupling of the two flavors. The previous analysis shows that the 
two contributions to Hh from the points go an d qi in momentum space decouple. In position 
space this tells us that a^ x splits up in the low-energy regime as a^ x = {X/i p |<i/e e tpx a qo+p }e zq ° x + 
{X)ipi<i/e e %px Q'q 1 +p}e iqiX ■ We write this instead as a^ x = a>(j t o)^ x e mx + a(j t i) tX e mx , in which the 
operators a(j^), x are given by the expressions in the curly brackets of the preceding equation. 

In the continuum limit we then have that e~ 1 / 2 a(j iAt ) iX tends to a smooth field ^f(j tll )(x), aris- 
ing from envelopes of plane waves around the point q^. This real-space description has been 
useful recently to generate non-trivial field potentials from lattice models [I]. For example, 
one might perturb Hh through the addition of an on-site potential such as e fi k (x)ol x a kiX . 
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k (x)e ^ 9 " )x a( jtlJ ,) tX a(k,v),x, however, if the function f jk (x) 
does not vary rapidly from site to site, the presence of the highly oscillatory term e~ l ^ qti ~ q ^ x 
will ensure that only fi = v will contribute in the continuum limit, and so the two flavors 
will decouple. At the other extreme, one can consider functions on the lattice that vary 
sufficiently rapidly, as say f(x) = g(x(qo — qi)), which can be used to produce non-trivial 
couplings of the flavors [JJ ITU], however here we do not consider such scenarios. 

3.4. Generation of the full transfer operator. It is straightforward to produce a mass 
term in the continuum limit, simply by adding the term J2 x m ( a },x a j,x ~ &LA>)> wm ch then 
generates the usual dispersion relation E 2 = p 2 + m 2 as e — > 0. However, we can more 
generally use functions {m J k (x, t)} that do not vary too rapidly over the lattice, and perturb 
Hh by the on-site potential term H m (t) = J2 X m f( x i t)( a ], x a k,x — b^ x bk,x)- This can be written 
as 

(3.7) 

H m{t) = { m o k ( x ^)( a lj^),x a (k^),x - h \j^). x h {k^),x) + e ±l(go ~ gi)x (flavor coupling terms) ) 

x,j,k,fi 

where the flavor coupling terms do not contribute in the continuum limit, as explained in 
the previous section. In addition to the terms Hh and H m for the auxiliary system we add a 
coupling term between the auxiliary and physical systems, which will generate the state |% e ). 
For this, we mirror the local coupling in [2] and define the coupling term 

(3.8) H mt (t) = R jk (^ t) [(«L a M + b lxh,x) ® fyx, y = t) 

X 

where again, for simplicity, we assume the functions {R^ k (x,t)} vary sufficiently slowly over 
the lattice so that the fj, = 0,1 flavors decouple once more. 

The final transfer operator that generates the PEPS state is finally given by 

(3.9) M e = 1 + e(H m + H h + H mt ) , 

where the operator hat serves to specify the non-trivial action on the physical field system. 
The basic tensor structure of M e {t) is shown in Fig. [TJ where there is also an implicit physical 
index at each site, coupling to the physical field A, which we omit in the diagram for the sake 
of clarity. 

In the next section we analyse the continuum limit, and derive the desired path integral 
representation from the smooth fields e~ x l 2 a^ ^ jX and e~ 1 ^ 2 by^^ x . Assuming that m J Q k and 

R? k slowly vary on the lattice, the expressions for H m and H mt in terms of these smooth fields 
are obtained by the doubling of flavor index j — > (j,fi). The same happens for the indices 
of the kinetic term Hh, but with added feature of a parity flip a x <H- <J z a x a z relating the two 
flavors. This discrete transformation is itself a Lorentz symmetry and ensures that the total 
field state will possess the full symmetry group, as we shall see later. 

3.5. Construction of the path integral representation. In this, and the subsequent, 
subsection we construct the path integral for the continuum limit of the sequence \Xe), £ 

0. We follow the one-dimensional prescription, and insert coherent-state resolutions of the 
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identity into the product TJ^ =0 [M e (t)]. This is first computed for a simple elementary 
timestep e and the continuum limit taken in the spatial direction. Finally the continuum 
limit is taken in the timelike direction, to obtain the final path integral representation of the 
cMPS state. 

For clarity we shall write to denote the compound index (j, /i) with j = 1, . . . , D and 
ji = 1,...,D for the two flavor sectors. This notation is useful since the actual PEPS 
parameters m?Q and R? k , that we can control do not have any \i dependence, and so the ji 
label simply doubles up the auxiliary fields, without playing an independent variational role. 
In what follows use coherent-state resolutions of the identity which, in both the fermionic 
and bosonic cases, are given by 

(3-10) \<f>{x,t)) = | ® w ,s0 Jm ,s(M)> := D({4> jli!S (x,t)})\n B ), 

where we have the usual coherent state displacement operator, given for a single mode with 
annihilation operator c as D(a) = exp[ac^ — ca*], the label s = a,b labels the particle 
type, and in the fermionic case 4>j s {x, t) are Grassmann numbers. We also use bold-faced 
4>(x,t) := {(f)j )S (x, t)} to suppress indices when the contractions are clear. The identity 
contribution is easily calculated via the overlap formulae for coherent states (see e.g. [2] for 
details) and gives (<f>(x, t + e)\<ft(x' , t)) = exp[—j^2 x (<p J! (x,t)d t <p(x,t) — d t (f>(x, t)^(j){x, t))], 
however the and H m terms require more attention. 

It is simplest to work in momentum space, for which 
(3.11) 

where m J k (p,t) = ^^ x e~ ipx m J k (x,t) anc ^ the indices on m^ fc and J jfc are explicitly inde- 
pendent of /i. 

Instead of expanding in terms of fermionic/bosonic coherent states of a,j x and bj )X we 
expand in terms of fermionic/bosonic coherent states of a,j p and bj p . Insertion of the above 
H h + H m into (4>(p,t)\H h + H m \4>(p',t)) gives, to 0(e), 

(3.12) n- 1 ' 2 ^ o fe (p-p^t)[0* MIa (p,^)^,a(p^^)-C, fc (p,^)^,6(p^^))+ 

s/N8(p - p')J^ k p'(^ a (p, t)4>^, b {p', t) + fj^ b {p, t)^ a { P \ t)] 
= N- 1 / 2 < k (P - *A t)<r z * k ,{p', t) + VN6{p - p')J^ k (t)(^l(p, t)p'a x ^^{p', t)) 

= £ mi k (x,t)^l(x,t)a z ^^(x,t) + J^t)^l(x,t)id x ^a x a^(x,t), 

where tyj^x, t) = ^^."'"[j^) J is a two-component auxiliary field with flavour index j^, and 
where j M = 1,2, ... ,D for each of the two separate \i = 0, 1 sectors. 
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The interaction term can be evaluated in the same way, and we obtain 

((J>(x,t + e)\H int \<t>(x',t)) = (^x,t)^l(x,t)^(x,t)^ A (x,y = t), 

where the hat on ip A is again to emphasize that the term is an operator on the physical 
system, as opposed to tyj (x,t) which is an auxiliary two component (Grassmann) spinor, 
with flavor index j M . 

The total sum in the spatial direction can now be evaluated, and becomes an integral over 
x, which provides us the single time-step amplitude coming from M e (t). Once done we can 
then sum over the time direction to obtain the final expression for the field state. However 
the process requires care, and in particular a field rescaling to preserve finiteness, and so we 
discuss this in the next section. 



3.6. Field rescaling and the continuum limit. For the graph used in the previous section, 
and also the square-lattice model in the appendix, the two-dimensional contraction across the 
graph requires a passage to the continuum in two independent directions, and so must be 
handled carefully. In this section we briefly spell out the technical details showing that we 
obtain a well-defined two-dimensional action, and we temporarily suppress the flavor-doubling 
label n to reduce our index load. 

Recall the basic form of the 2D PEPS state: 

lx)=C[A-... fci) ...A-.. feM) ...]|r 1 ...r M ...) 

= c[---]\n A ), 

where [%m • • • ku) are a set of contraction indices dependent on the particular graph structure 
of the state, and \r M ) is basis state at lattice site M. As we have already explained, the con- 
traction can be rewritten as a dynamical process involving the product of transfer operators 
generating infinitesimal steps 

(3.14) £[■■■] = (u L \M e (l)M t (l -€)■■■ M e (e)M e (0)\u R ). 

For clarity, we restrict to a finite set of points {(x, £)} where x runs over N x points, separated 
by a distance e x and t runs over N t points, in timesteps of e, and we make explicit all indices. 

Once we have introduced coherent state resolutions of the identity at each graph point we 
have at each value of t a total of 4DN X complex-valued functions to integrate over, coming 
from 

(3.15) 

k,s,x 

Inserting N t + 1 such resolutions into (I3.14p gives 

'l+e 

\(®k t+ ,, St+t ,x t+t (t>k t+ ^ t+ Xx t+e ,t + e)\M(t)\ ® kt 



(3.16) 



€[■■■] = J dfi 



i=0 



(uJ L \(l)k l + e ,s l+ X X l+eJ + e))(0fc O ,soOo,O)|WR), 
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with the measure for the integral given by 

, l+e 

( 3 - 1? ) d » = ^DN x (N t+ l) n 



t=0 



M.St.Xt 



The amplitudes in the integrand of (13.161) have been calculated in the low-energy sector and 
in this regime we obtain the expression 

dfJt e £^e[-±(*i M (x,t)9 t * fcM (x,t)-^ 

for the tensor contraction £[••-], where the matrices m and R are allowed to vary smoothly 
in both x and t. 

We now take the spatial e x — > limit followed by the temporal e — > limit. However, to 
keep things well-behaved, we must first rescale the integration variables 

(3.18) * kli (x,t) ^ -$=V kll {x,t). 



Indeed, this rescaling was to be expected in order to respect the correct commutation/anti- 
commutation relations: in 1+1 dimensions we have *f>k ~ a k / 'y/e^, while in 2+1 dimensions 
we should instead have ~ a^/ \Ze x e y , where e x and e y are the two spatial lattice scales. 

Once this rescaling is performed, we find in the e x , e — > continuum limit that the resultant 
cMPS field state becomes 

and we obtain a well-defined two-dimensional action. 

3.7. The two-dimensional variational class of field states. Once the rescaling has been 
performed we may directly integrate by parts, to obtain the class of two-dimensional field 

states 
(3.19) 

\ X (J,m ,R)) = J x?v^pvl,ty^(^ 



where we have introduced the coherent field state \<f>(x,y)) = exp[J dxdt(<&(x,t)ip A (x,y 

> ] - ( 



t) - $(x,t)ip A (x,y = t))]\£l A ) with $(x,t) = R jk (x,t)^/\ (x, t)^f k (x, t), and there is an 



implicit sum over /x, j and k. 

We have obtained this state through the continuum limit of a PEPS contraction, and so 
(13.191) describes a class of field states that inherits the desirable properties of PEPS, such as 
entropy area laws. The general action given above is not necessarily invariant under 5*0(1, 1). 
However, Lorentz symmetry for the field can be achieved for a subclass of states in which 
Ji k {t) = iS^ k . This allows the momentum terms to respect the desired symmetry and we 
obtain 



(3.20) \x(m,R)) = J P^p^ e i/^(^( i5 ^ t)), 
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where we have also let = —im? k . For the situation where m is diagonal and constant 
over the auxiliary spacetime, the coherent field amplitudes are then recognised as Grass- 
mann/complex number path integrals for a set of 2D uncoupled Lorentz invariant spinor 
fields. The spinors have the associated gamma matrices 7 = o z and 7 1 = ia y obeying the 
Clifford algebra relations {7^,7^} = 2if v \ with rf" = diag(l, — 1). Also for ji = 1 flavor 
sector we also perform a parity transformation on — > 7°^^, which is given by 7 = o~ z 
for the 1+1 dimensional case. 

Note that, in obtaining this symmetrical state, we were forced to take both m and J to 
be purely imaginary, which means the term Hh + H m corresponds to a generator of a unitary 
transformation on the auxiliary system. In particular U(l,0) = (1 — ieH(l))(l — ieH{l — 
€))■■■ — Texp[— i f H(t)] generates the cMPS state, where H(t) is built out of second-order 
and fourth-order combinations of creation and annihilation operators for the auxiliary and 
system fields. Furthermore, the distinctively two-dimensional term is the contribution from 
J(t), which couples the two spinor degrees of freedom. By setting J(t) = we reduce to a 
diagonal scenario, of the same form as obtained for the one-dimensional cMPS. 

Of course j M = vf^vl/ is a conserved current for the free Dirac field, and in particular 
is its charge density. Thus, in the case where the auxiliary state has manifest Lorentz 
symmetry, and where the different flavors decouple, we have the appealing interpretation that 
the matrix R, which is allowed to vary in both x and t, couples the densities of the different 
flavors and dynamically generates the physical field state. 

Of course, we could now weaken the conditions and allow more general J(x,t) and m(x,t) 
to obtain a Dirac action on a 1+1 dimensional spacetime with non-trivial metric. To do 
so in general would additionally require the modification of the identity term 1 to include 
a well-behaved function T(x,t). Looking back at the analysis, the key feature involved in 
the derivation is that the auxiliary particles have two degrees of internal freedom and are 
allowed to hop left or right with some amplitude or remain stationary while flipping an 
internal (spin) freedom. This is reminiscent of Feynman's 'checkerboard derivation' [HI EZ] 
of the Dirac propagator in 1+1 dimensions from a discrete lattice model. There, an electron 
moves along infinitesimal lightlike trajectories, while jointly flipping direction and spin under 
a Poisson process with an imaginary rate 1/im. In light of this, it is not so surprising that we 
have obtained a Dirac-like action in our continuum limit, although for us a key component is 
that an expansion in terms of coherent states in the auxiliary time direction either side of an 
operator M e = 1 + eH to t effectively generates a Legendre transformation of a 'Hamiltonian' 
H tot [7r, (/)]. For a coherent state-expansion, we obtain a term i<j)W t 4> from the coherent state 
overlaps which can be viewed as ird t (f) where 7r = i(jy is the momentum conjugate to 0. 

3.8. Analytic continuation to the euclidean sector. We now have a manifestly Lorentz- 
invariant auxiliary action, however the resultant physical state for A will have non-trivial 
entanglement structure in general. As such we would like to analytically continue to the eu- 
clidean sector and arrive at a rotation invariant auxiliary action. For second-order actions we 
can achieve this analytic continuation simply via t — > it, but for spinors subtleties arise. This 
coordinate transformation, when carried over to the Lorentz transformation, does provide the 
correct rotation group, but when acting on the spinors themselves results in no longer 
transforming as a scalar. A direct euclideanizing of fermion fields results in a number of 
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problems, such as a loss of hermiticity within the euclidean propagator. However, these diffi- 
culties were overcome by Osterwalder and Schrader [15] by making use of a construction that 
involves fermion doubling where the number of degrees of freedom are doubled so that the 
spinor and conjugate spinor are independent of each other and transform appropriately under 
the Lorentz group. However, there exist alternative approaches to euclideanizing the field 
that do not require this. Instead of analytically continuing the coordinates it is possible to 
analytically continue the metric itself if* — > r]^ v {9) so that it forms a one-complex-parameter 
family of metrics interpolating between the Minkowski and the euclidean one [IT]. A more 
abstract formulation can be achieved by using vielbeins, but for our purpose we work directly 
with the spacetime metric. 

The appropriate metric is rj^ y {9) = (cos 20/1 cos26 l |, — 1) defined for < 9 < n/2 ex- 
cept for the singularity at 9 = it /4, which can be circumvented via extending 9 to be 
complex. To carry the symmetry group over we also demand that the Clifford Algebra 
relation {7^(6*), jv(9)} = 2rj fJLV {9)l holds for all 9, and introduce the 75(6*) matrix obeying 
15{^V = 75(0)> {75(^)57^(^)1 = an d l5{9) 2 = 1. A particular parameterized set of gamma- 
matrices [H] that allow the continuation are then given by 

7i (0) = 7i 

(3.21) 70 W = |cosW /2 (7 ° CQSg + ' 75Sing) 

= 1 LlTrt(75Cos#-i7osin#), 

I cos 29\ 1 ' z 

where 75 = a x for our 1+1 dimensional case. 

This family of gamma matrices obeys the correct anti-commutation relations, and more 
importantly, generates an interpolation from the 5*0(1, 1) Lorentz transformations to the 
SO (2) rotations via A (a;; 9) = exp[u) tJ ' u 'E flu (9)], where oj parametrizes the group and 12^(9) = 
lv {9)] are its generators. The intuition behind this choice of parameterization is that 
to construct scalars under the Lorentz symmetry 5*0(1, 1), fermions in the (0, |) representa- 
tion are contracted with ones in (|,0), whereas for 50(2) we form contractions within the 
same representation [IS]. Since 75 is reducible over the different helicities, while 70 is not, 
the intuition is to rotate 70 -H- 75 to pass from 50(1, 1) to 50(2). 

A one-parameter family of actions, symmetric under this group action, can then be con- 
structed, where 9 = is the Lorentz-invariant Dirac action and 9 = n/2 is the desired 
rotation-invariant euclidean action. It is found to take the form 

(3.22) 5[^, f ; 6] = J rf 2 x v / -det(r / (^))$ t 7 O ^(^)7 At (0)^ - m]*. 

By inspection, we can see that this action could be obtained as the continuum limit of a 
one-parameter family of discrete tensor networks described by M e (t; 9) of the form 

- — -. . cos# sin# V cos 29 , ^ . 

(3.23) M e (t; 9) = t~ r = ~ <r y -j== + e A ^f ( H ™ + H h + H- mt ). 

Vcos2# vcos2# y | cos 29\ 

The parametrization of the metric has been transferred to the tensor contraction across 
the discrete graph, and we can smoothly transform from 6' = 0to6' = 7r/2to obtain, in the 
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continuum limit, the rotation-invariant euclidean action, and the cMPS class 
(3.24) \x{m,R)) = J exp[-Ss[W, * f ]]|$), 

where Se = f d 2 x^ {i^){i$ E — m)^ is the rotation- invariant action with 7^ = —270, and 
the flavor indices are implicit, and are the same as in equation (13.201) . 

The above analysis works independent of whether we have used bosonic fields or fermionic 
fields, however in the latter case it is possible to adopt a related approach that does not 
require the modification of the gamma matrices. Instead the Grassmann spinor and its 
conjugate \j/(x) are considered independent of each other, and instead of working with the 
full metric rf v one can analytically continue vielbeins e™ —> e™(9) = e l9 8™ and use these to 
construct a Dirac action in the way one would for curved spacetimes [26] . 

3.9. Area law properties. The two-dimensional field state \x) given by (I3.24p naturally 
inherits local properties from the discrete state. In particular, the resultant state is necessarily 
local in its entanglement structure and obeys an area law. For example, we could consider 
a finite region A of points on the graph described above. For this region we can define its 
boundary dA as the set of points in A within a distance e of points not in A, and \dA\ as 
the number of points in this set. For any pure quantum state \i/jab) there exists a unique 
measure of entanglement, namely the von Neumann entropy of entanglement, or simply the 
entanglement entropy, defined as Sa = — Tr[p^ logp^] where pA is the reduced state on 
system A. 

For discrete states the dimension D of the tensor labels within the contraction places an 
upper bound on the rank of the reduced state on A. The local definition of the state means 
that only systems in OA are entangled with the region B and so the entanglement entropy is 
upper bounded as Sa < c\dA\ where c is a constant dependent on D. For the continuum limit 
we need only impose that the region A is of fixed area with boundary of fixed length \dA\. 
Since the number of points in dA will scale linearly in 1/e, we then deduce that Sa < c\dA\/e 
and so it is clear that the resultant field state also obeys an area law. 

4. Conclusions 

In this paper we have constructed an abstract class of physically natural higher-dimensional 
quantum field states emerging from an underlying discrete PEPS state. Central to the con- 
struction was the recent path-integral representation introduced for one-dimensional cMPS 
states in [2] . The decomposition of the one-dimensional states into a superposition of coher- 
ent field states, together with an action for the auxiliary system, provided a useful guide in 
extending the class of field states to higher dimensions. 

The desired symmetries of the physical states are encoded in the auxiliary dynamics, and 
thus to obtain rotation-invariant states we may considered the imaginary time evolution of 
a Lorentz-invariant auxiliary field theory in (1 + 1) dimensions. To obtain this from first 
principles we began with a discrete PEPS on a particular graph, and demonstrated that the 
low-energy sector gave rise to a class of states which manifestly included rotation-invariant 
states. The continuum limit of the lattice state was shown to produce a doubling of flavors for 
which the two sectors can be made to simply decouple for PEPS data that varies sufficiently 
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slowly with respect to the scale of the network. Significantly, this doubling only occurs for 
the auxiliary degrees of freedom while the physical fields are left unaltered. The rotation 
invariance of the resultant field state arises from the symmetries of an auxiliary spinor action 
realised as the analytic continuation of a Lorentz-invariant action to imaginary time. Since 
PEPS states automatically obey area laws, we deduce that so too do the constructed quantum 
field states. 

There are no restrictions on the statistics of the physical fields, and since variational meth- 
ods evade the sign problem afflicting Monte Carlo techniques, it would be of interest to see 
how the cMPS variational class performs in finite fermion density scenarios. 

Finally, while in principle there is no obstacle to obtaining a continuum limit from a PEPS 
class with a bosonic auxiliary action with the desired symmetries, we leave this avenue to 
future work. 
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Appendix A. Continuum limit of the square lattice PEPS 

It is useful to spell out what exactly happens when we try to proceed directly from the 
square lattice PEPS to the continuum, and show that the resultant state is not rotationally 
symmetric, but retains features of the underlying lattice. 

We work with a square lattice of points labeled as (x, y) with the sides of the square of 
length e. To each vertex there is associated a tensor AL^-, with r the physical index for 
the resultant state and the remaining four indices denoting the two directions entering and 
leaving the point. The indices k, I take on discrete values 1, . . .D. 

We shall find that it is convenient to make use of the diagonal coordinates u — -(a; — y) 
and v = \{x + y) to describe the points in the plane. In these coordinates we have points near 
(x,y). Note that basic translations in this coordinate system are of the form (±e/2,±e/2). 
Here v is the auxiliary time direction, while u is the auxiliary spatial direction. 

The basic form of the cMPS state is 



where C[-} denotes contraction over all indices in accordance with the graph. 

We quantize the indices k G {1, . . . D} and (s = a,b), by mapping them to orthonormal 
states \k,s), and take these index states as orthogonal one particle modes of a bosonic sys- 
tem \k,s) — > \lk,s) '■= a !j^B)- F° r clarity, we can then view these index states as living 
in a Hilbert space H a (n,w) (g) "H fe (n,t;), associated to the lattice point (u,v). We consider 
v =constant lines, and define 



(A. 2) Uauxiv) := <g> ue zH a (u, v) ® Hb(u, v) 



where ultimately we will identify Uauxiv) and l-iaux{w) for any two times v and w. Also, for 
simplicity we consider a fixed array of indices B = {r(u, v), • ■ ■ ,r(p,q), ■ ■ • } defined over the 
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2D square lattice of points. Graphically, the contraction term is represented as 
(A.3) 



n\ A ri . . . A 



r-N 1 




A A 



A A 



A A 



A A 



A A 



where a contraction cut has been made along the v = constant line. Lines with arrows 
correspond to uncontracted indices, with right/up-pointing ones corresponding to bra- indices 
and left / down-pointing ones corresponding to ket-indices. The state of the auxiliary system 
at any fixed timestep is taken along similar diagonal lines, with boundary states \ojl) and 
\u)r) as shown. 

A.l. The contraction as an inner product. For any fixed diagonal v = vq we can split 
the tensor contraction £[•••] into two parts and write the result as an inner product £[•■■] = 
(u l (vo)\uj r (vo))(B), where (oj L (vo)\ denotes a 'forward pointing' tensor and \oj r (vo)) denotes 
a 'backward pointing' tensor (see figure |A~3|) . Furthermore, we have that 

\uj r (v )} G H aux (v ) 
(A.4) (u L (v )\eK ux (v ) 

We capture the dependence of £[•■-] on A^^ by writing \uj R (v )) = U[r(v )]\ui R (vi)) . 
Here we have advanced the cut from v = v q to v = v q + e/2 and have the new 'frontier' vector 
\co R (vi)) G WauxiyQ + e/2) and a transfer/transition operator 

(A.5) U[r(y )] : U aux (v Q + e/2) ^ H aux (v ), 

being a linear operator that describes the advancing contraction cut (or equivalently the 
infinitesimal evolution of the state \uj r (vi))), and r(v ) = (ri(v ), r 2 (fo), . . .) denotes a list 
of physical indices along the v — v o timeslice. 

The operator f/[r(uo)] is a string of tensor terms built from an operator M (different from 
the M defined in the paper). In particular 



(A.6) 



U[r(v )] = [M [n(u )] ® M[r 2 (v )] ® M[r 3 (v Q )] ■ ■ ■ ] 
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where crucially, M[r u (vo)} is a mapping between the Hilbert spaces 



(A.7) M[r u (v Q )} : U a (u + e/2, v + e/2) <g> - e/2, t; + e/2) -)• 7f a (u, u ) ® «o) 



and is given in terms of the PEPS tensors as 



(A.8) M[r u (v Q )\ = 1 + A^\ij)(kl\. 




where v n = v o + ne/2. 

A. 2. Construction of the Path Integral. We assume that at each point (u, v) on the 
square lattice we have the index state space T-L a (u,v) (g) T-L b (u,v), and the resolution of the 
identity operator for this space given by 



where \4>k, s ) = ex P[0M^l s — <fi*k s ^fc,s]l^s) = -D(0m)I^-b)- We also have the identity operator 
for %aux{v) as = ®«1( m ,d)- Here, we do not consider a single tensor product of Hilbert 
spaces at every point on the lattice, but instead view things 'dynamically' as a tensor product 
defined along a line v = constant, and so operators like U(r(v)) and l v are linear mappings 
from "H aux (» to H awc {v + e/2). 

To construct the path integral, we begin with 



and as before, we insert a resolution of the identity either side of C/[r(u)]. That is, we consider 
I( v+e /2)U[r(v)]l( v y This provides us with amplitudes 

(«><&,„(«, v) | <8> (®i^, 6 (M, u) |M[r u (u)] I ® fc fe , a (u + e/2, u + e/2)) | ®, ^ b {u - e/2, u + e/2)) . 
The central amplitude is then 
(®fc0fc,aO,^)| ®k <Pk,a( u + e / 2 ' v + e / 2 ))(®j0i,fe( M ^)l ®j <A?,&( M - e/ 2 ,^ + e/2)). 

+ A SfeSi)<®*^,a( u > v) | li,a> <®j^i,6(«, v) ] l ii6 > <l fc>a | <gi fc ,a(« + e/2,v + e/2)><l J>6 | <8> <j> ltb {u - e/2,v + e/2)) 
From the formula for the overlap of coherent states, the identity term gives the amplitude 



exp[-e/4(0 fcia (<9 u + d v )(p* k a + 4> k , b (-d u + d v )<f>% fi - (<f>* kia (d u + d v )(p k<a + <f>* kib (-d u + d v )<j> h)b ))] 
+A^ fe ( J o <0i ;O (u, v)|l ij0 ><^-, 6 (i*, ^)]l i)6 ><l fej£l |^ feja (u + e/2,v + e/2))(l,, 6 |0, j6 (« - e/2, v + e/2)) 



(A.IO) 




fe,s 



(A.ll) C[..-] 



(u; L ( W o)|t/[r"(wo)]^[r"K)]^[r( y2 )]t/[r( U2 )]---?7[rK)]|^K)) 



which can be written more compactly by defining ^/ k '■= (<j>k,a, (j>k,b), to give 



exp[~((^(a z d u + ld v )% - {^ k )\a z d u + ld v )V h )\ + 
+A \ik^ji)( ( l ) iA u ^ v )\h,a)(4>j,b( u i v)\lj, b )(l k , a \4> k ,a(u + e/2, u + e/2)) (1 I)6 | ^ )6 (u 



e/2, v + e/2)) 
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Note that transforming to the original coordinates, x = u + v and y — v — u we then have 
the first term contributing the amplitude 

(A.12) exp[i6S[y,**}} = exp[~(tfJV • % - • **)] 

defining an action S in either the (u, v) coordinates or in the original (x, y) coordinate system. 
It is clear that this action is not rotationally invariant. 

The above term is for the identity contribution in U. A particular tensor network state on 
the square lattice is described by and so by following the above prescription, we can 

create the appropriate operators within the one-particle sector of the Fock space. 

A simple and natural candidate is to take M of the form M = 1+ 1 (Q l J (u, v)a\bj <S> I& + I a <S> 
Q l £ ' (u,v)b\bj), where for example a>j,bk are single particle annihilation operators for mode-j 
in H a and mode-fc in H^, and {Q%(u, is a set of functions defined over the lattice 

points. The interpretation of this Q-term is that it generates scattering of particles within 
the internal bond degrees of freedom, without any scattering in the spatial degrees of freedom. 

Then for example, if we assume Q % 1 = Q(u,v) we automatically obtain the term 

(A.13) e -Q(u,v) + M) = \Q{u,v)^ 

i 

in the auxiliary system action after the above expansion in terms of coherent field states. 

Of course other options exist for the local contractions. For example, we can consider 
mixing the spatial modes with a term such as 

M(w, v) = 1 + eR IJ a\aj 
(A.14) = l + eR^^\4 a ^bl)(a Ja ®b lb ) 

where R IJ (u, v) can be taken to be a creation operator for particles in the physical field state. 
Clearly the addition of such a term generates quartic terms in the action of (<fr a , fa). 

The above analysis provides us with an amplitude contribution along a spatial slice, however 
to then sum over the auxiliary time and obtain the full amplitude requires a re-scaling of the 
field. The result of this is to yield an action 

(A.15) S = J dxdyh^\V ■ ^ k - Q(x, y)&V). 

We have shown the explicit details of this rescaling in section ( 13. 6ft . but there for the case of 
the rotationally invariant action. 
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